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Abstract 

We describe an application of computational techniques to the study of instanton induced 
baryon decay in high energy electroweak interactions. 
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We wish to illustrate an application of 
numerical methods to the semiclassical study 
of instanton induced baryon number viola- 
tion. In this investigation, computational 
tools take over where more traditional an- 
alytic methods reach their limit, and the 
two approaches complement one another. 
Our presentation will be in the form of a 
progress report since the computation offers 
formidable challenges and we have not yet 
brought it to completion. We believe, how- 
ever, that we have made enough progress 
to warrant the presentation of our interim 
results, both because of the intrinsic inter- 
est of the technique and in the hope that it 
may serve as inspiration for similar compu- 
tational investigations. 

Since the pioneering work of 't Hooft 
it has been known that electroweak inter- 
actions can give origin to baryon number 
violation because of the anomalous diver- 
gence of the baryon current and the occur- 
rence of topology changing instanton pro- 
cesses. Nevertheless, at low energy and low 
temperature the amplitude for baryon non- 
conserving processes is abysmally small be- 
cause of the barrier penetration factor: 

A ^ exp(-^i„,t//) , (1) 

where g is the electroweak coupling con- 
stant and Sinsti the action of the instan- 
ton configuration in rescaled units, is a nu- 
merical factor of order one. It has been 



shown that the barrier can be overcome by 
thermal fluctuations if the temperature be- 
comes comparable to the energy Esph of the 
sphaleron, the metastable state which sits 
on top of the barrier separating field con- 
figurations of different topology. A much 
more controversial question is whether un- 
suppressed baryon number violation can also 
occur in high energy collisions when E 
Egph- The difficulty stems from the fact 
that one needs to consider processes where 
the initial state is an exclusive two parti- 
cle state. So far semiclassical techniques 
have provided the only clues to the study 
of instanton induced baryon number viola- 
tion, but these methods are not amenable 
to the study of exclusive state processes. 
Recently a framework for solving the prob- 
lem has been proposed by Rubakov, Son 
and Tinyakov [§. (The literature on in- 
stanton induced baryon number violation is 
very rich, and although many authors de- 
serve credit, limitations of space force us 
to cite only work directly related to our in- 
vestigation.) The main idea is to consider 
processes where the initial state is also in- 
clusive, but constrained so that the particle 
number operator takes a fixed value = 
u/g^. The total energy is given hj E = e/g"^ 
and, with u and e held fixed, the small g 
limit allows one to evaluate functional inte- 
grals by a saddle point approximation. This 
leads to a cross section a ~ exp[— F(e, h')/g^], 
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and the the hmiting value of F for u 
provides information for processes with 
small numbers of initial particles. 

With some non-trivial manipulations one 
can determine the saddle point configura- 
tion that dominates the functional integral 
p|. This is given by the solution of the 
complexified classical equations of motion 
along a special contour in the complex time 
plane, shown in Fig. 1. In general such a 
solution cannot be found analytically, but 
its numerical calculation is within reach of 
today's computers. This is especially true if 
one considers either a simplified model, or 
configurations with rotational symmetry re- 
ducing the degrees of freedom to one-space 
and one-time dimensions. 
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Figure 1. Complex time contour. 



The boundary conditions are of special 
importance. The solution must be "real" 
along 7, where the asymptotic fields for t 
+00 represent the final state (if the system 
contains complex fields (p and (p, by real- 
ity we mean that = (p*). The condition 
on the particle number in the initial state 
translates into the requirement that for very 
early times in the evolution (i.e. asymptot- 
ically for t —oo + %T along a) the fields 
must reduce to a superposition of normal 
modes of oscillation with amplitudes satis- 
fying the equation 

Ai-k)=e'A*ik) , (2) 

where 6 plays the role of a chemical poten- 
tial conjugate to the particle number in the 
initial state. It is clear from Eq. ^ that for 
9^0 the fields cannot be real along a. 

Computationally, it is convenient to think 
of a time reversed evolution by which an 
initial configuration at t = undergoes a 
Euclidean evolution to t = iT along the 
imaginary time axis (i.e. following the ori- 
ented segment — and then a Minkowski 
evolution (along —a) from iT to — oo -|- iT. 
At t = 0, for an x-axis discretized with 
sites, we have independent real field coor- 
dinates and A^ pure-imaginary momentum 
components. The conditions on the normal 
mode amplitudes given by Eq. ^ amount to 
A^ complex constraint equations, or 2A^ real 
equations. Thus, in principle, one could 
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evolve the fields from an initial Ansatz at 
t = and adjust the initial variables so 
that Eq. || is satisfied. In practice, since 
the evolution equations along the imaginary 
time axis are elliptic, one cannot perform a 
forward integration. Rather, one must re- 
sort to some relaxation procedure or other 
global algorithm by which one solves the 
evolution equations set of simultane- 
ous non-linear equations for all points of a 
space-time grid. The situation is further 
complicated by the fact that with complex- 
ified fields one cannot just minimize a Eu- 
clidean action integral. We have developed 
a "second order" formulation, by which we 
minimize a constraint functional obtained 
from the modulus squared of the functions 
that must vanish at all grid points (an ear- 
lier part of this study was done in collabora- 
tion with Timothy Vaughan). The formal- 
ism of lattice gauge theory has been used 
to obtain a gauge invariant discretization. 
We tested our procedure in the context of 
the 2D Abelian Higgs model (one space and 
one time dimensions), where we found that 
it did reproduce the expected Euclidean so- 
lutions, including solutions with multiple 
bounces of the fields between two different 
topological sectors. 

Another crucial ingredient of the calcu- 
lation consists in extracting the amplitudes 
of the normal modes of oscillations once the 
fields reach the linear regime. As the energy 



disperses in the course of the Minkowski 
evolution, one expects that the system even- 
tually reaches a regime where the fields de- 
part little from their vacuum values. At 
this point the evolution becomes governed 
by the linearized field equations and quan- 
tities like normal mode amplitudes and par- 
ticle number can be meaningfully defined. 

Initially we studied the approach to the 
linear regime in the 2D Abelian Higgs sys- 
tem convenient simplified model which 
exhibits the anomaly and the non-trivial 
topological properties required for baryon 
number violation. We found, however, that 
with this system the decay of sphaleron like 
configurations gives origin to persistent, lo- 
calized excitations with an extremely slow 
rate of decay. These excitations are remi- 
niscent of the breather modes of 2D theories 
with real scalar fields that exhibit soliton 
behavior. The Abelian Higgs model has no 
soliton solutions in one spatial dimension, 
so one expects the localized excitations we 
observed to eventually dissipate. Still, their 
extremely long persistence made the com- 
putational extraction of the normal modes 
amplitudes quite problematic. 

This led us to consider the more realistic 
4D SU{2) Higgs system where the energy 
should disperse much more readily because 
of the higher dimensionality of space. We 
implemented the evolution equations for the 
spherically symmetric configurations stud- 
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ied by Ratra and Yaffe The spheri- 
cally symmetric configurations have a resid- 
ual Abelian gauge symmetry and are de- 
scribed by a two component U{1) gauge 
potential a^{r,t), and two complex Higgs- 
like fields x(^) ^) ^"^^ 4'{^jt) of charge 1 and 
1/2 respectively (x parameterizes the trans- 
verse components of the original 4D gauge 
potential and the original Higgs field). 
The decay of the sphaleron in this system 
has been studied by Hellmund and Kripf- 
ganz who observed the onset of the linear 
regime. 




Figure 2. Sphaleron decay in the four di- 
mensional SU{2) Higgs model: evolution of 
the X field. The values of the phase of the 
complex field are coded by different shades 
of grey. When the evolution reaches the 
linear regime, a gauge transformation, indi- 



cated by the sudden change of shading, is 
performed to extract the normal mode am- 
plitudes. 

We have developed code that simulates 
the Minkowski evolution of the above sys- 
tem in the ao = gauge. Fig. 2 illustrates 
the decay of a sphaleron after a small initial 
perturbation and the eventual onset of the 
linear regime. We have calculated the nor- 
mal modes both numerically, as eigenmodes 
of the discretized system, and by solving the 
continuum equations of motion analytically. 
There are 4 independent set of modes, cor- 
responding to 3 massive gauge bosons and 
1 massive scalar boson. In Fig. 3 we dis- 
play the time evolution of the total particle 
numbers of the 4 sets of modes. The con- 
stancy of these observables at late time is a 
good indication that the system has indeed 
reached the linear regime. 
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time contour is no longer required. Work on 
this problem is also in progress. 
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Figure 3. Sphaleron decay in the four di- 
mensional SU (2) Higgs model: behavior of 
the particle number in the four normal modes 
of oscillation of the linearized system as func- 
tion of time. 



The identification of the normal-mode am- 
plitudes also opens the possibility of study- 
ing the induced violation of baryon num- 
ber in the classically allowed region with 
energy above the sphaleron barrier. In this 
region transitions that change the topology 
are not suppressed by a barrier penetra- 
tion factor, but the question remains open 
as to whether they can occur in processes 
with an exclusive two-particle initial state. 
Semiclassical methods can be applied again 
(cf. Ref. 0) and the analysis is somewhat 
simpler because the detour along a complex 
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